Morphology of epitaxial core-shell nanowires 
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We analyze the morphological stability against azimuthal, axial, and general helical perturba- 
tions for epitaxial core-shell nanowires in the growth regimes limited by either surface diffusion 
or evaporation-condensation surface kinetics. For both regimes, we find that geometric parame- 
ters (i.e., core radius and shell thickness) play a central role in determining whether the nanowire 
remains cylindrical or its shell breaks up into epitaxial islands similar to those observed during 
Stranski-Krastanow growth in thin epilayers. The combination of small cores and rapid growth of 
the shell emerge as key ingredients for stable shell growth. Our results provide an explanation for 
the different core-shell morphologies reported in the Si-Ge system experimentally, and also identify 
a growth-induced intrinsic mechanism for the formation of helical nanowires. 



The combination of dimensionally confined trans- 
port and engineered interfaces inherent in heterostruc- 
tured nanowires has led to their emergence as an 
important cla ss of low dimensional, multifunctional 
nanostructuresfn^ In particular, radially heterostruc- 
tured core-shell nanowires (CSNWs) have enjoyed con- 
siderable attention from the synthesis and design commu- 
nities, mainly due to an unprecedented range of reported 
electronic properties such as high mobility carrier trans- 
port, tunable band ga ps, non-linea r optical gains, and gi- 
ant magnetoresistance! 4 * 5 * 6 * 7 * 8 * 9 * 10 * Modifications to usual 
nanowire synthesis routes, based on tailoring the rate 
and chemistry of radial growth during catalysis and/or 
vapor transport, have allowed the realizati on of several 
nanowire systems, including Si/Gep^ 12 * 13 ^ GaN/GaPP 
GaN/InGaNpl GaAs-GaJni^Ppl and metal oxides 
such as Sn0 2 /In20 3 P 

The direct control over the size and composition of the 
core and shell during synthesis should translate to wires 
with tailored orientations, degree of interface mixing and 
overall surface structure. Surface morphology is crucial 
as it drastically modifies the eventual function, yet its 
control remains an open issue. Indeed, recent studies 
in epitaxial systems have shown that the wires do not 
always remain cylindrical but can develop nanoscale sur- 
face modulations that tend to self-organize in a manner 
akin to quantum dots on thin films. For example, con- 
sider the Si-Ge system. The experiments of Lauhon et al. 
revealed that for small core radii (r c ~ 10 — 15 nm), the 
shell grows uniformly in both Si/Ge and Ge/Si nanowires 
(see Fig. Ilk) PI For large core radii (r c > 50nm) , Pan et 



* Corresponding author, email: mupmanyu@mines.edu 



al. found that Ge shells grow non-uniformly on Si cores 
and eventually develop into well-defined Ge islands ™ as 
shown in Fig. [TJd: this observation strongly suggests that 
the growth occurs in the Stranski-Krastanow mode in 
which the shell first develops a wetting layer and sub- 
sequently forms islands that grow and coarsen. Sim- 
ilar islands morphologies have also been observed in 
PbSe/Pb^niand Sn0 2 /In 2 3 CSNWsP 

The reports of these intriguing morphologies lead to 
the question: what conditions determine the stability of 
the epitaxial shell? In addition to the classical inter- 
play between surface energy and mismatc h strain that 
determines the stability of epilayers j 18 l 19 l 2Q l the shell mor- 
phology now also depends on core size and curvature, its 
growth direction and crystallography. In this letter, we 
use an isotropic continuum approach to study the com- 
bined effect of mismatch strains and geometrical param- 
eters on the morphology of epitaxial CSNWs. While we 
concentrate the comparison with experiments mainly on 
the Si/Ge system, the analysis we carried out is general 
and provides a fundamental framework for predicting the 
influence of geometry on the morphology of any epitaxial 
core-shell nanowire. 

To illustrate the main ideas, consider the scenario 
where a thin shell of thickness t = r s —r c and with average 
mismatch stress a s breaks up into n stress-free cuboidal 
islands of volume v — axaxh and average surface energy 
7 S . The energy of nanowire segment of length I is 

SF = -^{rl-rl)l + n(±ha ls ), 

where the first term is the energy gained via elimina- 
tion of the elastic strain energy and the second term is 
the increase in surface area due to island formation. As- 
suming that the islands self-assemble on the core surface 
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(Si/Ge) systems with different core radii. We propose that 
the classical interplay between the mismatch strain at the 
core-shell interface and surface energy depends on both the 
core (r c ) and shell (r s ), such that small fluctuations of the 
compressively strained Ge shells during growth can either de- 
cay, rendering the nanowire cylindrical (a,c), or can grow and 
cause the shell to break up into epitaxial islands (b,d). [Per- 
missions to reproduce the viewgraphs in panels (a) and (b) 
have been requested from the respective authors and publish- 
ers, and are currently pending.] 

into a hexagonal array with an inter-island distance A, as 
shown schematically in Fig. [T]i, the core area occupied 
by each island is A ~ A 2 . The surface area of the core 
is conserved, i.e. nA — 2itr c l. Setting SF(X) = 0, we 
find that the preferred inter-island spacing A m decreases 
parabolically with shell thickness, 

= J > * « rc (1) 

Using values for the Si/Ge system (7^ w 0.5 Jm -2 , 
E w lGPa, a « 5 GPa) and the reported island sizes and 
spacing in the unstable wireJ^ (a « 50 nm, h « 30 nm, 
A m ~ 100 nm), we extrapolate the maximum thickness 
of a stable shell to be t < 1 nm. 

A major oversimplification in the above analysis is that 
the average stress is assumed to be constant while in re- 
ality it also depends on the nanowire geometry (radii r s 
and r c ) and material parameters (core and shell lattice 
parameters a s and a c , Young's moduli E s and E c ). The 



analysis also does not make any connection with depo- 
sition conditions, in particular the surface flux. In the 
remainder of this article, we present results of a more 
rigorous linear stability analysis to study the morpholog- 
ical evolution of a dislocation- free, strained shell. Small 
perturbations in the otherwise cylindrical shell and the 
associated change in elastic stress distributions are ana- 
lyzed for stability with respect to surface relaxation ki- 
netics. The elastic solutions are quite detailed and we 
only summarize the main results in the text; the detailed 
derivations will be presented elsewhere [HW, MU and 
CC, in preparation]. For clarity, the results are limited 
to Si-Ge nanowires although the qualitative trends apply 
to any epitaxial CSNW system. 

Since our focus is on the effect of geometry, we (i) ig- 
nore intermixing at the core-shell interface, (ii) approx- 
imate the core and shell as isotropic elastic media, and 
(iii) assume that no facets exist at the core-shell inter- 
face. Then, the logical reference for the linear stability 
analysis is the elastic state of cylindrical CSNWs deter- 
mined by the dilatational mismatch strain (a s — a c )/a c 
at the core-shell interface, r — r c . For simplicity, we 
ignore differences in Poisson's ratio v of the core and 
shell. The strain distribution in the core and shell follow 
from the displacements (per unit length), which depend 
on the geometry and the mismatch. Ignoring non-linear 
elastic effects, Hooke's law and eq uilibri um conditions 
readily yield the stress distribution! 21 * 22 ^ Axial symme- 
try dictates that the variations in stresses and strains 
are limited to in-plane normal components within the 
shell. Therefore, the surface stress state is entirely due 
to the non-zero normal components with similar magni- 
tudes, and (j* zz . Taken together, they set the energy 
per atom on the surface, £* . We must emphasize that 
this surface elastic energy decreases rapidly with the shell 
thickness as it approaches the dimensions of the core, a 
trend that we shall later show to have a dominant effect 
on the morphological stability of the shell. 

To make the stability analysis more transparent, we 
first consider azimuthal perturbations. The nanowire 
surface geometry now takes the form r{6) = r s +A cosn#, 
where is the (in-plane) angular coordinate, A is the 
perturbation amplitude, and n = 27r/(A/r s ) > 2 is the 
azimuthal wave number. Geometrical constraints limit 
the possible wavenumbers to n > 2. The solution to 
this plane strain problem is the unperturbed, cylindrical 
core-shell system superposed with appropriate boundary 
conditions on an interior shell surface defined by the pre- 
scribed perturbation. The stress state is extracted using 
Airy stress functions for the core and the shell (Mitchell 
solution^ an d imposing interface conditions at r = r c . 
The traction free conditions at the surface are absorbed 
by transforming the stress tensor into local normal (77)- 
tangent (r) frame 19 , and setting a vv = a vr = 0. The 
equations are solved for small slopes and amplitudes, 
nA/r s < 1. 

The resultant stress driven morphological evolution is 
analyzed at the shell surface with the volume diffusion ig- 



3 




FIG. 2: (a-d) Dispersion relations uj(n) for the azimuthal perturbations of the shell when the growth kinetics is limited by 
surface diffusion (black curves), or by evaporation-condensation (gray curves). The growth rates uj are in units of (Bt)- 1 . The 
geometrically forbidden regime of n < 2 is shown as a vertical gray stripe in each case. The dispersion relations for the Si/Ge 
systems with t — 1 nm thick shell are shown in panels (a) and (b) for core radii of r c = 50 nm and r c — 10 nm, respectively; 
panels (c) and (d) show the dispersion relations for Ge/Si systems with the same geometrical parameters as in (a) and (b). 
The inset in (b) illustrates the stabilizing effect of increasing the shell thickness by a factor of 10. 



nored. The central quantity which drives the relaxation 
is the surface chemical potential, which we define with 
respect to that of the cylindrical nanowire, /i*. In or- 
der to capture the essential physics, we limit ourselves to 
an elastocapillary potential, ji s = ji* + rtj s tz + (£ — £*), 
where Q is the atomic volume in the shell, 7 S is the sur- 
face energy, n is the surface curvature, and 8 is the misfit 
strain induced elastic energy (per atom) on the surfaced 
While non-linear effects induced by wetting potentials 
and surface stresses can certainly become important pi 
they are the focus of follow-up work and ignored in this 
study. We limit ourselves to the two limiting cases of sur- 
face kinetics, i.e. surface diffusion and vapor transport 
via evaporation-condensation at the surfaceP^^ In the 
absence of an external deposition flux or at small depo- 
sition rates, surface diffusion is expected to control the 
relaxation (SD-limited kinetics). The flux in this case is 
proportional to gradient in chemical potential at the sur- 
face, V sH s . In the presence of an external surface flux, 
as in shell growth via chemical vapor deposition (CVD), 
the morphological evolution is limited by evaporation- 
condensation at the surface (EC-limited kinetics). At 
near equilibrium, the flux is proportional to the differ- 
ence between ambient and local surface vapor pressure, 
which in turn depends on the local curvature. 

In both scenarios, the surface stress state of the per- 
turbed nanowire modifies the surface flux and therefore 



the radial growth of the shell. In the small slope limit, 
the governing partial differential equations yield a pertur- 
bation amplitude which varies exponentially with time, 
A(t) = A(0) exp[Bcut], where B is a positive constant 
that depends on the material properties, deposition con- 
ditions and surface kinetics.^ Clearly, it is the sign of oj 
which determines whether the perturbation grows (un- 
stable kinetics) or decays away (stable kinetics). For the 
azimuthal perturbations, this reduced growth can be ex- 
pressed as 

^ = H^h^ + ^^ M T (2) 

where M is a positive function that depends on n, ge- 
ometry, and material properties. Since the general func- 
tion is rather complicated involving powers of n, in the 
remainder of this study we investigate its effect numeri- 
cally for the Si-Ge system. The function f(n/r s ) = /(A) 
is positive and depends on the type of surface kinetics; it 
scales as A -3 and A -1 for SD- and EC-limited kinetics, 
respectively. The remaining expression (inside brackets 
of Eq. [2| is the critical condition since its sign determines 
the stability of the system. The first term captures the 
stabilization due to high surface energies and surface cur- 
vature (k oc l/r s ), while the effect of elastic stresses is 
embodied in the second term. Setting uo = 0, we arrive 
at the value of the critical perturbation wavenumber n cr , 
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which is identical for the two classes of surface kinetics.^ 
A direct comparison, however, cannot be made between 
the two as the growth of the perturbation amplitude A(t) 
also depends on deposition parameters, lumped in the 
constant B, which are different for the two cases ™ 

Quantitative insight into the nanowire stability can 
be obtained from the variation in the reduced growth 
rate with the perturbation wavenumber, uo vs n, or the 
dispersion relations. The plots are shown in Fig. [2] for 
Si/Ge and Ge/Si systems and for two different core radii, 
r c = 50 nm and r c = 10 nm. The choice is based on 
the nanowire geometries reported in Refs. fTTI and fT2l (see 
Fig. [I]). In each case, the results are reported for a shell 
of thickness t = 1 nm. We find a dramatic effect of the 
core radius on the range of unstable wavenumbers. In 
the case of the Si/Ge system, the critical wavenumber 
for the larger cores is n cr « 23 (Fig. [2^l) . The corre- 
sponding azimuthal diffusion length is A « 14 nm. The 
prefactor /(A) leads to surface kinetics dependent maxi- 
mal growth modes n m (value of n for which duo/dn = 0, 
also the most unstable mode), quite like the trends re- 
ported for stressed thin films ! 19 * 20 ! For EC-limited kinet- 
ics, the maximal wavenumber is n m = 10, and the az- 
imuthal (diffusion) wavelength associated with this insta- 
bility A m ^31 nm. This length is smaller for SD-limited 
kinetics, X rn « 17 nm. 

Figure [2}d reveals that the range of unstable wavenum- 
bers and the corresponding reduced growth rates uo are 
significantly smaller than those corresponding to r c = 

50 nm (Fig. . While the diffusion lengths do not 
change, they are now comparable to the core radius, 
approaching half the circumference of the core cross- 
section. For example, n m w 3 during SD-limited kinet- 
ics (Fig. J2Jd), resulting in large relative diffusion length, 
A m /r s ~ 2. We expect the kinetics to be further re- 
tarded as the diffusion lengths are now spread over much 
larger variation in core curvature such that atomic-scale 
effects such as barriers at surface steps start to become 
important. 26 The same trends are predicted for the Ge/Si 
system (Figs. [2^ and[2]i), illustrating again that at con- 
stant shell thickness t a small value of core radius is a 
key factor for the stability of the cylindrical shell. The 
parameters considered in Fig. [2]i (r c = 10 nm, t = 1 nm) 
are quite similar to those corresponding to the Ge/Si 
CSNWs recently synthesized by Lu et a/P^, r c — 7.5 nm 
and t =2-5 nm. Using the experimental values, we have 
found that the diffusion wavelength is about two times 
larger than the core radius. Our analysis thus predicts 
that the Si shell of Ge/Si CSNWs is stable, which is con- 
sistent with the experiments of Lu et al 13 The effect of a 
stiffer Si shell under tensile strain is to increase the range 
of unstable wavenumbers and the corresponding growth 
rates. Again, however, strain induced thermodynamic 
and kinetic effects are not factored in our analysis and 
they could become important for stabilizing the tensile 

51 epilayers in these nanowires. 27 

A similar approach is employed to investigate shell 
stability with respect to axial perturbations. The shell 



surface morphology depends on the axial coordinate, z. 
We have considered perturbations of the form r(z) = 
r s + A cos kz, where k is the axial wavenumber. The 
elastic solution can be expressed in terms of Papkovich- 
Neuber potentials. 23 Since the core and shell have differ- 
ent elastic moduli, the relevant Airy stress function now 
involves modified Bessel functions of the first and second 
kind. The linear stability analysis reveals that the dis- 
persion relations, uo vs kr s: are similar qualitatively to 
those corresponding to azimuthal perturbations. This is 
not surprising and merely indicates that trends dictated 
by the balance between surface and elastic energies, cal- 
culated in the isotropic limit, are approximately the same 
for axial and azimuthal perturbations; the effect of inher- 
ent anisotropy in geometry on stability of thin shells that 
we have considered is relatively small. 

The inset in Fig. [2}d shows the dispersion relation for 
azimuthal perturbations in the Si/Ge nanowire for a ten- 
fold increase in shell thickness (t = 10 nm). The shell 
has the same dimensions as the core, t = r s . Clearly, the 
increase in shell thickness has a remarkable effect on the 
stability. The growth rate is now negative for all possible 
perturbations indicating that at least for the small core 
sizes employed in the experiments of Lauhon et al., there 
exists a critical shell thickness beyond which the shell 
becomes stable. In order to ascertain if this is a general 
trend, we have systematically analyzed the stability of 
Si/Ge nanowires with various core and shell sizes. The 
results are presented in the form of stability diagrams 
- contour plots of the critical azimuthal (n cr ) and axial 
wavenumbers (k cr r s ) as functions of nanowire geometry 
(Fig. [3]). For large core radii, the sharp turning point 
(knee) in the contours indicates that thin shells (r s /r c 
ratio below the knee) become progressively unstable ini- 
tially. Interestingly, further growth past the knee reverses 
this trend. The enhanced stability is more pronounced 
for small core radii, where the contours no longer exhibit 
a turning point. Increasing the shell thickness always in- 
creases its stability with respect to both azimuthal and 
axial perturbations. 

It is informative to discuss the predictions of stability 
phase diagrams such as those in Fig. [3] for other sys- 
tems than Si/Ge CSNWs. We have found that for the 
Sn0 2 /In 2 3 (Ref. 16) and GaN/GaP (Ref. 14) nanowire 
systems the experimental observations are in close agree- 
ment with our theoretical stability analysis. The phase 
diagrams for these two systems are similar to those given 
in Fig. H for Si/Ge nanowires. The GaN/GaP system 
has a uniform cylindrical stable shell: this is consistent 
with our theory, which based on the geometric param- 
eters of Ref. HH (r c = 7 nm, r s /r c = 3.0) predicts that 
the GaN/GaP system is deep inside the stable region 
for both azimuthal and axial perturbations. For the 
Sn02/In203 system}^ there are two distinct morpho- 
logical regimes observed experimentally, both of which 
agree with the theoretical analysis. The Sn02/In203 
CSNWs with smaller core radii and larger r s /r c ratios 
(e.g., r c = 33 nm, r s /r c = 1.76) are cylindrical and have a 
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stable shellp^ while the shell of the nanowires with larger 
core radii and smaller rJr c ratios (e.g. r c = 45 nm, 
r s/ r c — 1-64) is unstable^ and develops islands similar 
to those reported by Pan et aZP^ Lastly, we note that 
for Ge-core/Si-shell systems, the softer Ge cores result 
in Si shells which increase their stability irrespective of 
the core sizes considered (< 100 nm), i.e. there is no 
knee in the contours. The results of the stability anal- 
ysis are again consistent with the observations of stable 
shell Ge/Si CSNWs by Lu et a/. 13 , since the experimen- 
tal parameters place these Ge/Si nanowires close to the 
stability boundary (not shown) in our calculations. 

The trend is opposite to that predicted by Eq. [I] 
(A oc 1/vT), and also to that normally observed in com- 
pressively strained thin films on planar substrates! 19 * 20 * 
In these core-shell heterostructures with finite-size core, 
the elastic stresses and surface (shell) curvature both de- 
crease rapidly with increasing shell thickness. The stabil- 
ity diagrams show that the former is more sensitive as the 
shell approaches the dimensions of the core. The destabi- 
lization due to decrease in surface curvature (~ 7 s /r s , the 
first term in the parenthesis in Eq. |2| is relatively small, 
and overall the shell becomes more stable. The behavior 
underscores the importance of a finite core in absorb- 
ing the mismatch stresses in these low dimensional het- 
erostructures, thereby dramatically changing the balance 
between misfit and surface stresses that is at the heart of 
the Asaro-Tiller-Grinfeld instability in thin filmsP^ 

The stability diagrams provide an elegant yet sim- 
ple explanation for the variations in the morphologies 
of the as-synthesized CSNWs. To directly make contact 
with the Si-Ge system, we have overlaid the approximate 
growth path of the Ge shell in the two studies shown 
in Fig. [I] For the small Si core employed by Lauhon et 
a/.f^the window over which the compressively strained 
Ge shell is morphologically unstable is small (below the 
hashed region) and limited to early stages of growth, 
which explains the uniform, cylindrical shell morphology 
observed in these CSNWs. In addition to the balance 
between surface curvature and the elastic stresses, both 
of which decay with shell growth (Eq. [2j, we expect sta- 
bilization due to the kinetic barriers to the surface re- 
laxation that now needs to occur over larger lengths and 
surface curvatures. On the other hand, CSNWs with 
thicker Si cores reported by Pan et a/P^ result in a larger 
window of shell thicknesses over which it can develop 
axial or azimuthal instabilities. The range of unstable 
wavenumbers during early stage growth correspond to 
smaller diffusion lengths. The morphological instabil- 
ity is further aided by the reduced changes in curvature 
spanned by the smaller diffusion lengths, and also by the 
higher temperatures (~ 770 °K) employed for the shell 
growth. The latter follows from the fact that the uncat- 
alyzed CVD-based synthesis route for Ge shell deposition 
exploits the lower thermal stability of the precusor gas 
(GeHU), as opposed to the lower downstream tempera- 
tures (~ 700 °K) prevalent during stable shell growth in 
the study by Lauhon et al. The combination of these 
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FIG. 3: Contour plots of (a) the critical azimuthal wavenum- 
bers rtcr and (b) the critical axial wavenumbers k cr r s associ- 
ated with surface diffusion-limited relaxation of a sinusoidally 
perturbed Ge shell in a Si/Ge nanowire, as functions of the 
core and shell radii, r s and r c respectively. The hashed region 
indicates a stable shell. Note that the definition of a stability 
with respect to axial perturbations is somewhat arbitrary as 
the perturbation is limited by the length of the wire. The 
cut-off wavenumber below which we consider the nanowire to 
be stable corresponds to a nanowire length of approximately 
1 /im. The vertical dotted lines correspond to growth on Si 
cores with radii r c reported in Refs. [TT1and[T2l as indicated. 



factors makes the shell susceptible to island growth dur- 
ing early stages of its growth, which is precisely what is 
observed in the experiments.^ 

Closer examination of the reported island morpholo- 
gies in the unstable Si/Ge nanowires reveals an axial 
inter-island distance in the range A m = 75 — 150 nm. 
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FIG. 4: Contour plots of the reduced growth rate u(n,kr s ) a function of the azimuthal and axial wavenumbers for Si/Ge 
nanowires with two different geometries: (a) r c — 50 nm, r s — 51 nm, and (b) r c — 10 nm, r s — 11 nm. The growth rates are in 
given units of 10 -2 (Bt) -1 . The plots in (c) and (d) show the ratio a of the maximally growing helical to axial growth rates, 
a — u; m (l, kr s )/ujm(0, kr s ), for the k values corresponding to the points shown by the arrows in (a) and (b), respectively. The 
vertical dash lines in (c) and (d) correspond to the r s /r c values that were used in panels (a) and (b), respectively. 



The axial diffusion wavelength length associated with 
the most unstable mode predicted by our analysis is 
A m = 31 nm (EC-limited kinetics). The quantitative 
disparity is hardly surprising as we ignore effects such 
as intermixing and formation of misfit dislocations at 
the core-shell interface, growth crystallography, surface 
stresses, and faceting of the core. Faceting in particular 
appears to be quite important, as the islands are ob- 
served to form preferentially on planar facets; the stress 
relaxation at the corners results in enhanced stabiliza- 
tion, consistent with Eq. [2] 

We have also analyzed the interplay between the two 
perturbation modes, by considering general perturba- 
tions of the form r(#, z) = r s + Acos(n# + kz). Figure [4] 
shows combined effect of the two modes on SD-limited re- 
duced growth rate oj in the Si/Ge system, again for a shell 
of thickness t = 1 nm. The semi-circular contours follow 
from the fact that the effective wavenumber is a linear 
combination of the two modes, k e ff = \/k 2 r 2 + n 2 . For 
large cores (Fig. |4^l), the most unstable mode is clearly 
along the axial direction. Superposition of an azimuthal 
component always decreases a;, indicative of an intrin- 
sic stabilization mechanism that is likely to be activated 
due to the spatial variations that naturally exist during 
the uncatalyzed deposition of the shell. The extent of 
decrease depends on the magnitude of the scaling fac- 
tor for the growth rate, B. While the islands do form 
axially in the corresponding experiments,^ they cannot 
serve as validation since the Si core is distinctly faceted 
and therefore influences the interplay between the two 
perturbation modes. 

For smaller cores, the reduced growth rate exhibits a 
maximum for the general perturbation n = l,/cr s ~ 3 
(Fig. |4Jd). The helical radial growth of the shell, shown 



schematically in the figure, is simply a reflection of the 
energy balance captured in Eq. [2] The capillary and elas- 
tic energy terms are functions of the effective wavenum- 
ber such that, per unit length of the nanowire, addition 
of a helical component to an otherwise axial perturbation 
decreases the curvature dependent surface energy and/or 
increases the stored elastic energy in the shell. The ge- 
ometrical window associated with these unstable helical 
perturbations changes significantly with the growth of 
the shell, primarily due to the variation of the elastic en- 
ergy (term M in Eq. [2]) with the effective wavenumber 
k e ff. As an illustration, Figs. [4 z and|4]i show the effect of 
an azimuthal component (n = 1) on the maximally grow- 
ing axial perturbation function of the geometric 
ratio r s /r c . For larger cores (Fig. Efc), the ratio of the 
reduced growth rates a = uj m {l,k'r^)/(j m (0,kr 8 ) during 
initial stages of shell growth is slightly less then unity, 
consistent with the stability diagram reported for t = 
1 nm. As the shell thickness increases, we observe a win- 
dow of radii ratio for which helical perturbations should 
dominate the growth morphology, 2 < r s /r c < 2.7. In 
the case of the smaller cores, the ratio is greater than 
unity during early stages. However, the helical growth 
mode switches to a pure axial instability past a criti- 
cal radii ratio, r s /r c « 1.6: this variation underscores 
the importance of geometry on the evolution of the shell 
morphology. 

Krill et al. have predicted the stability of such he- 
lical morphologies via surface diffusion, albeit for pure 
whiskers subject to extrinsic strains. 28 In the case of epi- 
taxial core-shell architecture where the mismatch strains 
are present intrinsically, it may become prohibitive for 
the core to sustain the helical growth of the shell elas- 
tically. A likely scenario, especially when the shell is 
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softer than the core, is that the core buckles helically 
to better accommodate the perturbations of the shell. 
The net twist entailed by the buckled core can also be 
relaxed by defects such as screw dislocations. The fact 
that such morphologies can grow due to mismatch strains 
(and possibly surface stresses) appears to be a novel in- 
trinsic mechanism for stability of helical morphologies in 
core-shell nanowire systems. 

In summary, the results of the linear stability anal- 
yses for epitaxial core-shell nanowire systems highlight 
the importance of geometry in determining the nanowire 
morphology. The dependence on the core size allows us 
to rationalize the cylindrical (stable) and island (unsta- 
ble) growth observed experimentally in the Si-Ge sys- 
tem. The variation with shell size shows the importance 
of the growth rate, relative to surface relaxation kinetics, 
in controlling the shell morphology. Since the temper- 
ature and partial pressure of the precursor gas directly 
affect the growth rate, they become important deposition 
variables that allow control over the nanowire morphol- 



ogy. There are, for example, CSNWs systems 17 whose 
morphology can be changed from stable shells to epi- 
taxial islands solely by adjusting the growth rate, even 
in regimes of geometric parameters where our theory 
strictly predicts uniform shells. In future studies, we will 
envision the development of a more comprehensive un- 
derstanding of the morphological evolution of core-shell 
nanowire heterostructures in which the effects of inter- 
mixing, anisotropic surface energies, surface stresses, and 
kinetics could be addressed in various combinations. 
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